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AN INVESTIGATION OF THE EFFECT OF TWIST ON THE COUPLED 


BENDING-TORSION VIBRATIONS OF A CANTILEVER BEAM 


SUMMARY 

An experimental study was made of the effect of twist on the natural 
frequencies, coupled in bending and torsion, of a cantilever beam. The shape of 
the beam was selected to provide a high degree of coupling. An analytic procedure 
for determining these frequencies was also developed and is included in this report. 

It was determined that coupling without twist reduced all frequencies from 
their uncoupled values. The fundamental frequency was found to remain relatively 
constant as the beam was twisted. All natural frequencies higher than the 
fundamental were lowered as the total twist was increased to 15 degrees, but 


remained relatively constant with further twist. 





INTRODUCTION 

With the advance of helicopters and turbo-machinery, the problem of twist 
and its effect upon the natural vibrations of a beam becomes of more and more 
interest. Much work has been done to determine the effect of twist on beams 
having no coupling between torsion and bending, and the effect of a high degrec of 
bending-torsion coupling on untwisted beams. This report shows primarily the 
effect of twist on the natural frequencies of a particular beam highly coupled in 
bending and torsion. It also shows the node lines for all the natural frequencies, 
and indicates the effect of bending-torsion coupling on the natural frequencies in 
the untwisted beam. 

A beam with dimensions as shown in Fig. 1 is considered. This particular 
shape was chosen to give a maximum amount of bending-torsion coupling, 
relatively easy tooling problems, and readily obtainable sectional and elastic 
properties. In this report, the term coupling refers to the interaction of the 
torsional, flapwise bending, and chordwise bending types of vibrations. It was 
planned to use both experimental tests and analytic means to determine the 
natural frequencies resulting from the combination of twist and coupling. Only 
experimental results however, are contained in this report. Difficulties en- 
countered in preparing the digital computer program prevented the inclusion of 
analytic results. The mathematical method used is outlined in the Appendix. 

This method is based on one being used in a report of Isakson and Eisley 
currently being prepared for the University of Michigan Research Institute. It 


involves the introduction of coupling terms from the equations of motion developed 





by Houbolt and Brooks, Ref. 1, into the matrix method of solution of Targoff, Ref.2 &3. 
The Targoff method was described and used in a report of Isakson and Eisley, Ref. 4. 
The coordinate system used is shown in Fig, 2. The flexural center of the 
cross section was located according to a formula from Roark. Ref. 5, which formula 

was developed from Ref. 7. The elastic axis of the beam was maintained as a 
straight line by mechanically twisting the beam about this axis. The nomenclature 
and sign convention Íor the displacements, shears and bending moments used in the 
analysis are shown in Fig. 2. All parameters used or needed are listed in the table 
of symbols. The non-dimensional forms of these parameters used in the analysis 
are also listed in the table of symbols. The mass and stiffness distributions are 
uniform for the length of the beam. 

In gathering the experimental data, tests were made at total angles of twist 
BEN 5. 5°, 7.90, 12.99, 20.79, 26.79, 30.49, and 40,80, 

In conducting this investigation, the assistance and guidance of Professors 
Eisley and Isakson, of the Department of Aeronautical Engineering of the University 


of Michigan, were of great value. 
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Bending moment in flapwise direction 

Bending moment in chordwise direction 
Torque about elastic axis at any cross seotion 
Beam length 

Shear load in flapwise and chordwise direction 
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mode constant for determination of frequency 


distance between area centroid of tensile member 
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acceleration due to gravity 
polar radius of gyration of oross-sectional area, 
effective in carrying out tensile stresses, about 
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length of blade segment 
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total loadings per unit length in x N, and $ 
directions 

total torque loadings per unit length in x, N, an d f 
directions 
thickness of cross section at any chordwise station 
coordinate in direction of X axis measured from 
the root along the elastic axis 
blade angle between major axis of cross-section 
and the "y" axis for £? z O, 
positive when leading edge is up. 
rate of twist of blade per unit length 
increment in ß between blade segments 
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displacement of elastic axis in ''y' and " Z '' direction, 
respectively 
displacement of elastic axis in direction of minor and 
major principal axes of the cross-section, respectively 
total twist in blade between x = © and X = R 
cross-sectional coordinate; N axis lies along major axis; 
€ axis is perpendicular to major axis and passes through 
elastic axis. 
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angle of twisting deformation, positive when leading edge is up 


mass per unit length of blade 
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Note: all zero subscripts refer to root section 
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Parameter Values for Beam Under Consideration 


Values Independent of Twist 
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Non Dimensional Parameters 


Values independent of twist 
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Values dependent on twist 
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EQUIPMENT AND PROCEDURE 

The model used in the present investigation was made of 2024 T 4 Aluminum, 
machined from 6" x 1" stock. It was machined to within . 015 inch of the given 
beam dimensions. 

Four inches at the end of the model was left rectangular for mounting purposes 
(Fig. 1). 

The bar was mounted with four bolts between two modified channel sections, 
which in turn were fastened to a double box section bolted to a structural member 
in the wall of the building. The double box section was made by welding 1/2" inch 
steel plates between the outer edges of the flanges of a 6 x 6, 5/8 inch steel I beam. 
(Figs. 3 and 4). The channel sections were later reinforced by welding a 0.5 inch 
steel gusset as shown in Fig. 3. This was done to observe the effect upon the 
natural frequencies of substantially stiffening the mounting. No measurable change 
in frequencies was noted. This provided some indication that the mounting was 
sufficiently rigid to give results that closely approached a true cantilever. 

MB Vibration Test Equipment, Model T1 - 32034 was used to vibrate the 
beam. This machine has a frequency range of from two to 20,000 cycles per second. 
The vibrator itself was mounted to the floor and fastened to the beam with an 
aluminum tube. This tube was fastened to the vibrator with a flexure joint and to 
the beam with a universal joint (Figs. 3 and 4). 

Two universal joints were used originally, one at each end of the aluminum 
tube, but this arrangement allowed too much slack in the connection. Two flexure 


joints were also tried in the same manner as the universals, but this arrangement 





was too stiff to permit the beam sufficient freedom in torsion. 

To obtain the desired twist, special fittings were made to fasten to each end 
of the bar so that it could be twisted about its elastic axis. A Riehle torsion testing 
machine of 120,000 in. lbs. capacity was used to twist the beam. The twisting rig 
as shown in Figs. 5 and 6 was used to transmit the torque. 

The method of determining natural frequencies consisted of the following steps: 

1) A frequency was approximately located by listening for and observing the 
amplitude maximum. In the case of the fundamental frequency this method provided 
an accurate answer by noting the frequency at which the measured maximum amplitude 
of the free end of the beam occurred. 

2) The beam had been designed so that uncoupled natural frequencies of torsion 
and chordwise bending occurred in the vicinity of a natural uncoupled flapwise 
bending frequency. When the beam was being vibrated in the vicinity of a band of 
natural frequencies, fine sand of silicone flour was scattered over the surface of the 
beam. The frequency was then adjusted until sharp node lines appeared and concided 
with maximum amplitudes of vibration. The maximum amplitudes were recognized 
by first turning down the amplitude control of the vibrator until the sand or powder 
was barely agitated, then adjusting the frequency until the maximum agitation of the 
particles appeared. This method permitted the determination of amplitude maxima: о 
in a band of frequencies where sharp node lines did not disappear, but only shifted 
their positions. Where several maxima . appeared in a band, the strongest was 
assumed to be the frequency of the flapwise bending mode. This assumption was 
made considering that a flapwise node was designed to occur in each band, and that 


due to the shape of the cross section, the flapwise vibration would be the dominant one. 
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It is the frequency change due to coupling effects as influenced by twist 
rather than the absolute value of the frequencies that is important in this report. 
Several factors occur in the experimental procedure which would tend to produce 
an experimental result different from one arrived at by calculations. Throughout 
the experimental runs, attention was directed toward keeping these factors constant. 

The connecting rig between the beam and the shaker added end mass effect; 
this is discussed in the next section of this report. It was noted that varying the 
tightness of the bolts at the mounting even slightly could produce a slight variation 
in frequency. Of course, for all runs the beam was mounted as securely as possible 
so that the net effect of variations from one test to another was very small. 

For the evaluation fo the analytic part of this report the Royal Precision 
Electronic Computer LGP-30, manufactured by the Royal McBee Corporation was 
to be used. A matrix program was set up as shown in the Appendix. The computer 
evidently could not handle the problem as it was programmed. A great deal of 
time had been spent on this part of the investigation and insufficient time remained 


to pursue it further. 
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RESULTS AND DISCUSSION 

The effects of twist upon the natural freauencies of the beam are shown in 
Figs. 7 through 12. In Figs. 7 through 9, the absolute values of the natural 
frequencies are plotted versus the angle of twist. In Figs. 10 through 12, instead 
of absolute frequencies, the ratio of each frequency at any degree of twist to the 
corresponding value for the untwisted beam is plotted versus the angle of twist. 

The frequency ratios are grouped according to the assumed predominant motion, 
i.e. flapwise, chordwise or torsional. 

Table I lists the uncoupled natural frequencies of the beam calculated on the 
basis of geometry and mass of the beam. These are compared with the corresponding 
natural frequencies observed when vibrating the untwisted beam. A coupled “resonant 
frequency appears in the frequency spectrum corresponding to each calculated un- 
coupled natural frequency. It was assumed that the type of motion associated with 
each uncoupled natural frequency, remained as the predominant motion in the 
corresponding coupled natural frequency of the untwisted beam, and would continue 
to predominate as each frequency changed with twist. 

The discussion of the effect of twist is based primarily on Figs. 10 through 12, 
the frequency ratio plots. In this discussion each frequency is identified by its 
predominant motion. 

Twist has very little effect upon the fundamental frequency (first flapwise 
bending) of the beam. This frequency increases very slightly as the twist is increased 
through 30 degrees. The data point obtained at 40.4 degrees of twist indicates that 
there is a possibility of a more rapid rise in frequency as the twist is increased 
beyond this point. 
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The second flapwise frequency drops steadily as twist is increased to 
approximately 15 degrees; it remains constant as twist is increased further. 

The third and fourth flapwise frequencies decrease as twist is increased 
to 10 degrees. For the next 10 degrees of twist, these frequencies rise slightly, 
then remain constant as twist is increased further. 

The first torsional frequency, corresponding to the fourth natural coupled 
frequency of the beam, drops off rather rapidly as twist is increased to 15 degrees, 
it then decreases slowly as twist is further increased. The second torsional 
frequency decreases until about 15 degrees of twist is obtained, then tends to 
remain constant. 

The third torsional frequency deoreased to 15 degrees of twist then appeared 
to rise as twist was further increased. Data for the third torsional frequency was 
considered less reliable than the other data, because of the difficulty encountered 
in obtaining its resonance point and in obtaining sharp node lines. 

Only the first chordwise frequency, corresponding to the third natural 
coupled frequency of the beam, was obtained. It dropped sharply until a twist of 
about 20 degrees was obtained; it then appeared to remain constant or rise slightly 
as further twist was introduced. 

The shape and location of the node lines for the different frequencies are 
shown in Figs. 13 through 20. These were obtained when determining the natural 
frequencies. Although in some cases a band of frequencies occurred in which the 
node lines remained sharply defined and only changed their orientation and position 
with frequency, resonant peaks within the band were strong enough to permit 


determination of resonant frequencies. The occurence of frequency bands is 
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indicated in Figs. 13 through 20. 

The shape and location of the node lines was significant in that they helped 
in the identification of the type of motion. It was noted that the node lines tend to 
move toward the free end of the beam as twist is increased to about 20 degrees. 
With further twist they remain relatively constant or even move slightly away 
from the free end. The position of the node lines in sach group of frequencies 
reflects the flapwise bending mode occurring in that group. 

The node line illustrations show the appearance of what was assumed to be 
resonant frequencies of plate bending. This assumption was made since no 
natural frequencies based on flapwise bending, chordwise bending, or torsion, were 
predicted to oocur in this range. This plate bending first appeared near 300 cycles 
per second and became more pronounced as twist was increased. At high angles of 
twist, (20 to 30 degrees), it began at approximately 276 cycles and produced 
distinct node lines through a band of 400 oycles per second. Twist had little effect 
on the mean frequency of this mode as can be seen in Fig. 21. 

This plate bending appeared again at about 680 cycles for a twist of 25.7 degrees, 
becoming evry distinct at a twist of 30.4 degrees. Another plate bending mode was 
seen to oocur at 960 cycles and first appeared at a twist of 5.5 degrees. 

Fortunately, the plate bending resonant frequencies appeared between the 
other bands of frequencies and though they undoubtedly influenced the shape of the 
node lines in the upper range of frequencies, Å is felt that they did not greatly 
hinder the accurate determination of resonant frequencies in modes of vibration 
which are of interest in this paper. 


In shaking the beam, oertain frequencies above the fourth natural frequency 
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did not give distinct nodes. In these cases, the natural frequencies were plotted 

as a band rather than as a distinte point (Figs. 7 through 9). This difficulty was 

due not only to the proximity of the flapwise, chordwise, and torsional uncoupled 
frequencies, but also to the appearance of plate bending which was discussed above. 
Other inaccuracies were due to some unavoidable looseness in the vibrating 
mechanism, primarily the universal joint. This looseness also limited the frequency 
range obtainable, since the vibrating movement was lost in transmittal at high 
frequencies. The tube and fittings which connected the bean to the vibrator also 
added a slight end mass effect to the beam. By adding additional known mass 
weights to the end of the beam, it was determined that the effect of the vibrating rig 
reduced the fundamental frequency by 1.6 cycles per second. The data taken was 

not corrected for this effect, however, since the effect at higher frequencies was 
not known. The frequencies might also have been slightly low due to the impossibility 
of obtaining a perfectly rigid support at the root. The base of the beam however, 


was considered sufficiently rigid so that this effect was minor. 


15 





CONCLUSIONS 


E In all cases, the effect of bending-torsion coupling, independent of twist, 
appears to reduce the natural frequencies of the beam. 

2. The fundamental frequency is relatively insensitive to twist, demonstrating 
only a slight rise as the angle of twist is increased. 

3. All natural frequencies above the fundamental one are lowered as twist is 
increased to approximately 15 degrees. 

4, Beyond 15 degrees of twist the natural frequencies are relatively independent 


of further twist. 
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APPENDIX 
METHOD OF ANALYSIS 

The analysis is basically the one which is outlined in Ref. 3 and the 
Appendix of Ref. 4. It is extended to include torsional vibration and bending- 
torsion coupling and is applied to the non-rotating case. 

The beam is divided spanwise into ten equal segments. The mass of each 
segment is assumed concentrated at its center, and the bending stiffnesses EI 1 
and EI, and angle of incidence 8 are assumed constant between masses. The 
twist of the beam is accounted for by relative rotations of adjacent uniform bays 
(between masses) about a spanwise axis, the change in angle, A 8 , being equal to 
the total twist in a segment and occuring just outboard of the mass. 

With torsional motion included, the following quantities, at any point along 
the beam, may be expressed as a column matrix: 
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The elements of this matrix vary from one station to the next along the beam 
in the following fashion: LẠNG = IR | IE | IF | (A) „ where İRİ IE | [3 
are 10 x 10 matrices representing linear relationships, described below, between the 
elements of (Aİ and corresponding elements in (A) net Coupling effects are 
introduced into these relationships by consideration of the differential equations of 
coupled motion developed by Houbolt and Brooks in Ref. 1. 

In the [E] matrix which relates the { A) matrices across a vveightless 
section, the elements are obtained as follows: 
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For the untwisted section f > O , but the second term in the above equation 
is not dropped since it represents the coupling between chordwise bending 
and twist which can only appear at this piont in the relationships. For this 
term, P must equal the amount of twist that actually exists in a segment 
having length x. 

(E ese EB, gg)" = M” + VA x (1) 
(g) que = qu 


From Ref. 1, at any station: 
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4\2 
where A = GJ + EB 
NA 
Equations (a) through (j) provide the elements of the EJ matrix; these 
elements are put into non-dimensional form and the result is the 10x10 


| E] matrix whose terms are all zero except for the following: 
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The quantities B l and B 2 are those defined by Houbolt and Brooks, Ref. 1, 
and were evaluated for the particular beam. 
The elements of the [F] matrix which relate the (2) matrices on either 


side of a concentrated mass are determined as follows: 
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The second term is one resulting from the entire segment's undergoing a 
change of slope when the beam bends. 


k, = radius of gyration of a segment about a chordwise axis; m = mass of a 


segment. 
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The second term is a rotary inertia term as before; 


pe radius of gyration of a segment about a vertical axis; m=mass of a 


segment. 
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where Ko is the polar radius of gyration of 
the segment. 
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Equations (a) through (j) are put into non-dimensional form and provide the 


elements of the [F ] matrix: 
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The [R] matrix serves to rotate the coordinate axes through the angle A 0 


and appears as follows: 
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The analysis then continues in the same manner as in Ref. 4, except that the 
matrices are now 10x10 and yield a 10x5 product when the boundary conditions at 
the tip are introduced. 

The final results are obtained through evaluation of a 5x5 determinate. 


Introduction of trial values of N | D= f (w) ] , produce points 
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which define the residual curve. Zero values of this curve define the natural 


frequencies. 





TABLE I 
COMPARISON OF COUPLED AND UNCOUPLED FREQUENCIES 


I —— a 


Frequency Type of Vibration 
(Cycles per second) 
Uncoupled* Coupled 





lst Flapwise 


ist Torsion 
2nd Flapwise 
lst Chordwise 


2nd Torsion 
3rd Flapwise 


3rd Torsion 
4th Flapwise 


Does not include the effect of end masses. 
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Fig 1 


Sketch of Beam Showing Dimensions 
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